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ABGERACT 4; 

Underwater sound was scattered from a wind-agitated random 
rough model sea surface: The surface energy spectral density, 
the surface rms height, and the surface correlation function were 
measured. The mean square magnitude and probability distribution 
of the specularly scattered sound were measured for various angles 
of incidence between co and grazing. The results are compared 
with the basic theories expressed by Eckart and Beckmann and 
shadowing corrections derived by Wagner. The frequency spectrum 
of the scattered sound was studied and found to contain doppler 
shifted components, contrary to theory which predicts no doppler 
shifting in the purely specular direction. A detailed analysis 
of the origin of this doppler shifting is presented, and it is 
shown that interaction between axial and off-axis radiation is 


the explanation of the observed frequency broadening, 
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Specular Scatter of Underwater Sound 
from a Wind-Driven Model Sea Surface 
..” “ints Oductaen: 
Backgr ound 

When underwater sound is incident on the surface of the water, 
it is scattered. This scattering is dependent upon the character- 
istics of the surface at the time. If the surface is smooth, the 
scattering will be pure specular reflection; that is, angle of 
incidence equals angle of reflection; both are measured with re- 
spect to the normal to the smooth surface. If the surface 1s rough 
the scattering 1S no longer purely specular, but energy 1s Scat- 
tered in all directions. Two directions of scattering are parti- 
cularly interesting, back scatter and specular scatter. 

In recent years, since the mid 1950's, a great amount of work 
has been done in studying the back scattered energy. This case is 
relatively easy to study at sea, and many experiments have been 
performed. In 1967 Ball and Carlson studied this problem in the 
laboratory at the Naval Postgraduate School under controlled cir- 
cumstances. They used a water surface which was roughened by 
blowing air over it. Using this reproducible surface they examined 
the relative mean acoustic pressure and doppler shift of the back 
scattered energy. Their experimental apparatus will be referenced 
throughout this report. 

Distribution of the Scattered Sound Pressure 

Our experiments were conducted to measure the effects of water 

waves on scattering of underwater sound in the specular direction. 


The applicable theories on the subject are based on the Kirchhoff 


Es 


Assumption, which asserts that plane reflection coefficients can be 
used on a point by point basis for a gently curved surface. The 
basic theories assume that there is no shadowing or multiple scat- 
tering. Wagner (1966) has calculated a correction for the shadowing 
assuming ray propagation. This correction becomes important near 
grazing incidence. Our model experiment includes an oceanographic 
study of the scattering surface and permits a check of both the 
Kirchhoff-based theories and Wagner's correction. 

The experiments to measure the distribution of the scattered 
sound pressure were conducted by LCDR scheible and the results were 
presented at the May 1968 meeting of the Acoustical Society of 
America in Ottawa, 

Spectrum of the Scattered Sound 

Observations of scattered sound pressure have shown it to be 
both amplitude and frequency modulated. The most recent attempts 
to describe this behavior have been made by L. Wetzel (1966) and 
B.E. Parkins (1967). Although Wetzel's descriptions are explicitly 
derived for long EM waves, his conclusions are directly applicable 
to the acoustic problem. Wetzel bases his developments on a re- 
sonant interaction between the incident energy and an appropriate 
frequency component in the surface. Parkins uses physical optics 
approximations to describe the reflections occuring at the surface. 
Both of these developments arrive at predictions of doppler shifting 
of energy upon scattering from the surface. They predict doppler 
shifting for all directions of scatter except the specular. They 
thus imply that the energy scattered in the specular direction would 


be constant in amplitude and that it would not contain any frequencies 
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other than the incident frequency. However, it has been observed 
that specularly scattered energy does fade and may contain frequency 
or phase modulation. Mellen (1964) has observed this frequency be- 
havior for back scattering. The frequency spectra of specular 
scattered sound has been measured at sea, near Bermuda, by Roderick 
and Cron (1968). :° Their observations show the frequency spread to 
be less than 1 Hz under all ocean conditions and less than ,2 Hz 
under normal conditions. 

The experiments to measure the spectrum of the scattered sound 
pressure were conducted by LT Fowler and the results were presented 
at the May 1968 meeting of the Acoustical Society of America 


(Supplementary paper D3A). 
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2. Theory: 
Spectrum of the Scattered Sound 

The following treatment is a summary of Parkins (1967). "When 
a monochromatic acoustic wave is incident upon the time varying 
surface of the ocean, the envelope of the reradiation varies in 
time or fades." Parkins claims that this fading is a result of 
doppler shifting of the scattered pressure. His contention is 
that the incident radiation is "locally shifted in frequency" by 
scatterers moving on the surface. He assumes the surface to be a 
stationary, homogeneous, Gaussian process that is described by a 
three-dimensional spectrum. 


The locally incident pressure is given by 





-1ky “rae 
+; 1 1W t 
G9. =@ e 
i 
eA ; 
Ky = — (alee am + ‘kn) 


where 1, m, and n are direction cosines. Using developments 
similar to Beckmann (1963), the scattered pressure in the far field 


as a function of position and time is given by 


p 
° 1(Wt = kRo) q Hi aaae ‘ 
9 (p,t) = ~~ ~—__** ff eK * ax ay 
a -4 


27 R 
O 
where k = eo. 
iN 
A = acoustic wave length 
K=K,-« 


- k 
2 L 
area of integation = a square 21 ona side 


and r are defined on Figure 2.1. 
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To find the spectral density of the scattered pressure, Parkins 
transforms its auto-correlation function 
[ee] 


S(Q) = if R(t yeu ar 


nee) 





x 
where R(T) = (f(p,t)B (p,t + 7T)), that is the ensemble average, 


which is the time average if you assume ergodicity. He then shows 


: 2 
that for a slightly rough surface, K, o “<< | 


ie 2 52 mee Ree 2 BD 
S(Q)= B exp(-K_ 0 rE (Q-wW)16 Sin Kf Sin K ble K,, 
Zee 22 
Se) 
+ MOTKS O° G(K, 4K, ,0 ){ 
+ 2o> = 2 
where B = (@ /2nR ) (K-k, /K_) 
is 1(Q Ww )T 
G(K_,K_ ,Q-wW) = { CA @ Gary 6 tyer! ) dt 
x’ y oe x 
© ©) 
t = -7) 
G Wasi ES FO) ff p exp [ i(K u +; ey) du dv 
OX 
u = x!-x 
a ie 
p = surface auto-correlation function as a function of space 


and time 


O 


surface rms height 
Assuming a Neumann-Pierson directional wave spectrum and 
integrating, he then shows the spectral density of the scattered 


sound pressure as 
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$(2) = B exp(-k*o*)[ (2) 16 sinks Sin*k 8 /K IK 
Deo 2 2 > Bs 
+, 19% Kio SAN K 5K, )1 6 Ov - (KL 4K) 9“) 


+[ At( “Ky 5 -K,)1° 6 Q “w+ gKs) 49°) ] 


The importance of this equation is that is shows the reradiated 
pressure to be composed of a specularly reflected component at 
the same frequency as the incident pressure and two components 


that are doppler shifted by 


1 


2 4g? 


(Ko + Ke 

Parkins says, "The scattered components are perhaps sur - 
prisingly, given by delta functions."' He further states, "It 
seems likely that the spectrum of the scattered radiation will 
be narrow but not have delta-function behavior. The spectrum of 
the specularly reflected component depends only upon the gravity 
Wave spectrum through the phase modulation term exp ( -K* ), and 
it will consist of a single line at 2 =W as long as the sea, as 
seen from the angles of incidence and observation, is only slightly 
rough."' By slightly rough he means the roughness factor Ke Be 9g 
of Beckmann is less than 1. 

It should be noted that Parkins obtains both positive and 
negative doppler shifts by considering scattering from opposite 
sides of the origin, that is +K., ta: Mellen (1964) measured 
the spectrum in the back scatter case and found doppler shifting 
in only one direction. Further, his doppler shift is about one- 


half that predicted by Parkins. 


i 


In 1966, Wetzel, by considering a resonant interaction between 
incident HF radiation and surface wave lengths, also predicted 
doppler shifts. His scattering geometry is shown in Figure 2.2. 


He developes the relationships 


= + in 9. - Sin Q 
Ky tke St ; Sin 0. Cos oe 
= + in Q 1 
K, ye oa Sata! 3 Sin 2.) 
2 Zz 2 
K = Ky + Ky 
where Ky = the x component of the wave propagation number 
K, = the y component of the wave propagation number 
K = the wave propagation number 
OA 
XA = wave length of the radiation 


He argues that for conical scatter, Qo. = OL; or — 2. <1, the doppler 


shift 1s given by 


L. 1. 
- 2 4 
2, = (gk Sin ©,)* [2(1 - Cos g.)] 


where g = the acceleration of gravity. This equation also predicts 
zero doppler shift in the specular direction. 

Both of these developments have been for plane waves, that is, 
beam pattern divergence has not been considered. The assumption 
of a plane wave front is justified when the surface is far from the 
source and only a small portion of the wave front is considered. 
However, it seems likely that this will not always be true when 


considering sonar energy interacting with the surface. 
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Distribution of Scattered Sound Pressure at 
Near Grazing Incidence 
The magnitude of the specularly scattered sound pressure is 
a function of the relative roughness of the surface. The most 
commonly used measure of roughness was defined by Beckmann (1963) 
as the "roughness factor'' g, which is defined by the equation 


1 
2 


g°“ = ko (Cos eo, + Cos 9 where k = 27/A is the propagation con- 


2) 
stant of the sound, 9 is the rms wave height, and oat and 5 are 
the angles of incidence and reflection. For specular scattering 
Zhe a, ; 
oe = K. @9 ain Parkin's notation. 
An equation for the relative mean square pressure of specularly 
scattered sound, measured relative to mirror reflection, has been 


developed by Medwin (1967): 


ToL? A Cos 9 


(BB *) =e 9 + S(g). 
Yo 64 xo A* 


In that paper S(g) was defined, and it was shown to be approximately 
equal to g for g < 0.3. For our experiment, the second term, which 
represents incoherent scatter, is always less than 2% of the first. 


Therefore it has been neglected, and we have used e 9 


as the pre- 
dicted value of relative mean square pressure. 
Wagner (1966) derived a shadowing function making use of ray 


geometry and statistical mathematics. For specular scattering this 


function 


Ze 





ama v = Cot 0/2(0/L). 

Oo = rms wave height 

L = surface correlation distance 

Q = angle of incidence or reflection 
Wagner also shows in that paper that the correct application of 
this correction is to multiply the predicted value of relative 
mean square pressure by the square of this shadowing function. 

Medwin following Beckmann also developed an equation for the 
variance of the magnitude of the relative pressure for specular 


scattering: 


_ Tig a LA 
8 f\9, | - — 
2s” X xy 


Substituting the definition of g above and c/f for ik gives 


Sudt 2 16 L 

_ 27°O L Af Cos 0 

of ial ~~ 2 
1 


From this equation it can be seen that the variance should depend 
on the sixth power of the frequency and on the fourth power of 


the cosine of the angle of incidence. 
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3. The Model Sea Surface: 
General 

The physical arrangement of the experiment is shown in Figure 
3.1. A large anechoic water tank was covered with 4'* plywood 3%" 
above the water surface. Air, from Gee centrifugal fans, was 
forced between this covering and the water to generate a model sea 
surface. At the opposite end of the tank we have an "tanechoic" 
beach made out of six layers of small-mesh screen. 

In order to compare the behavior of sound pressure that is 
scattered from the surface with theory, one must know the surface 
energy spectrum, the surface correlation function R(d), the sur- 
face rms height 0, and the surface rms slope ¥. 

These developments usually assume a NeumannrPierson surface 
spectrum and a correlation function af either or 7 (a/b)? 
where the correlation distance L is defined as the distance at 
which the correlation is reduced to one e-th of its value at zero 
separation. 

For the surface used in this research the surface energy 
spectralfdensity, the spatial correlation function, and the surface 
rms height were measured. The frequency spectrum and rms height 
were measured at several positions on the surface (Figure 3.2). 
The rms slope was previously measured by Ball and Carlson. The 
rms height and slope were calculated from the energy spectral 
density and compared with measured values. 

The experiments that were conducted to measure the surface 
frequency spectrum, rms height, and correlation function were per- 


formed jointly by LCDR Scheible and LT Fowler, 
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Energy Spectral Density 

Electronics. Figure 3.3 shows the circuit used to measure 
the energy spectral density E(f) of the surface. The oscillator 
has a stable 10 kHz, 1.5 volt output. A resistive probe (Figure 
3.5) is placed in series sista the output impedance of the oscil- 
lator. This probe was designed and analyzed by Ball and Carlson. 
As the wetted length of the ah ole changes, its resistance varies 
and produces an amplitude modulation of the oscillator output. 
This modulation is then an electrical analog of the water surface. 
It is necessary to use at least a 1 volt oscillator output to 
ensure the modulation will be in the linear portion of the de- 
tector (Figure 3.4). Keller (1968) also used this detector and 
found it to be linear when the input is in excess of .5 volts. 

The detector output is the electrical analog of the surface. 
This analog has about a 1 volt DC component and a 10 kHz ripple. 
The differential amplifier (linear down to DC) is used to remove 
the DC and ripple. The DC is se by subtracting an equal 
value, oe set on the variable DC supply. The differential 
amplifier has an internal, variable, low-pass filter which is set 
at 200 Hz. The filter is flat from DC and down 3 db at 200 Hz. 
Because the surface frequencies are less than 100 Hz, they are 
passed by the filter but the 10 kHz ripple is blocked. The dif- 
ferential amplifier output is then recorded at .375 ips. The tape 
recorder is used in the FM mode to be sure the frequencies from 
DC to 100 Hz are preserved. The —— of the recorder is DC 
ver 10 Hz + 1 db, S/N 35°db. The tape recorder 1S increased in 


speed by 10, or 100 for play back to increase the surface 
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Figure 3.4 Detector 





Figure 3.5 Surface Probe for Measuring Wave Heights 
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frequencies to values that can be analyzed. This speed-up is 
necessary because the surface spectrum is predominantly from 0 to 
20 Hz, and the wave analyzer can not accurately measure frequencies 
below 10 Hz. 

The wave analyzer is used to measure the spectrum of the tape 
recorder output. When the tape recorder is played at 3.75 ips, 
the 10 Hz filter is used. In the first experiments performed, a 
complete spectrum was obtained for both replay speeds. [It was 
found that the predominant surface energy was at 3.8 Hz and that 
one need only analyze with a replay speed of 3.75 ips. This allows 
accurate spectral analysis of surface frequencies in excess of 1 Hz. 

The output from the wave analyzer is a time averaged voltage. 
With the analyzer set to a given frequency, it was found that the 
energy in that component varied with time. This variation was 
partially removed by using a slow response on the analyzer and 
was eliminated by eye averaging. The analyzer voltage was con- 
verted to a 1 Hz bandwidth by multiplying by the appropriate fac- 
tor (265 for 3 Hz BW, .375 for 10 Hz BW, and .159 for 50 Hz BW). 
The analyzer voltage is converted to ae by multiplying by the 
system calibration. The square of this result is the energy den- 
sity at that frequency with a 1 Hz BW (cm~-sec). A third correction 
factor 1S necessary to account for the multiplication of frequency. 
When the surface analog is replayed at 3.75 ips and the wave an- 
alyzer is set for a 3 Hz BW, the real frequency BW is in fact 
-3 Hz. Therefore the wave analyzer output energy, (voltage X BW 
correction X conversion factor)”, must be multiplied by 10. Like- 
wise, when the replay speed is 37.5 ips, the energy must be mul- 


tiplied by 100. 
30 


Calibration. To calibrate the output voltage as a function 
of wave height, the probe depth was changed in increments and the 
DC output from the tape recorder noted. To verify the linearity 
of the tape recorder, the DC voltage input to the recorder was 
also recorded. Figure 3.6 shows a sample of these results. A 
calibration was obtained for each spectrum obtained. Whe conver - 


wt 


sion factor 1S equal to 


A voltage DC - 1__inches _ A volts Zms 
A probe height (inches) 2.54 cm cm 


Although the voltages recorded were DC, the displacements were 
also DC; therefore, the results can be expressed as volts rms to 
cm rms. Ball and Carlson verified a static calibration to be 
valid for frequencies up to at least 100 Hz. 

Results. The surface energy spectral density was :determined 
for stations l C; 1 R, and 2 C of Pigure 3.2. These a are 
presented in Figures 3.7 through 3.9. Spectra were determined with 
and without the beach installed. It was found that the spectra 
without Pe beach were higher and peaked to the left of the ones 
with the beach. This is the result of reflection from the end of 
the tank. Since the objective of this report is to compare the 
scattered sound pressure with a sea-like surface, it was desirable 
to have the surface propagation vectors in the downwind direction 
only; therefore, it was decided to use the beach. 

Initially the spectrum was obtained with the probe attached 
to the plywood tank covering. It was felt that this plywood might 
vibrate and introduce noise into the system. Although this was 
not found to be true, the probe mounting was shifted to the metal 


grating which covers the tank. 
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Figure 3.6 Surface Probe Calibration 
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The system noise level was found to be down 10 db at 1 Hz, 
40 db at 4 Hz and 20 db at 10 Hz. The spectrum of the noise was 
determined using the same procedure as the surface spectrum except 
that the fans were off. 

Examination of the spectra in Figures 3.7 through 3.9 yields 
the following. The energy at station 2 is higher than at station 
1 by about 1 db, which is within experimental accuracy. The energy 
1.5 feet off the center of the tank is down by about 5 db, which 
indicates that the air from the fans was funneled to the center of 
the channel. The effect of the beach was to remove the energy 
reflecting from the end of the tank, thus lowering the total energy. 


All spectra with the beach installed peak at 3.8 Hz, have an g77°3 


-5.5 


dependence near the peak reducing to f 


ee 


creased frequency. The f 


dependence with in- 
dependence is in good agreement with 
the Neumann-Pierson spectrum and Phillips (1966) theory which 
predicts eo 
Using the surface energy spectral density it is possible to 
find the rms height 0 and the rms slope %&. These two parameters 


are used in various theories to predict the behavior of the scat- 


tered sound pressure. 


Integration of Figure 3.8 yields ao of 0.22 cm, which is in good 
agreement with that found using a sampling technique for the same 
surface as shown later in this section. O.M. Phillips (1966) 


shows that Bes 
oe af K~ B(w)dw , 
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io) 


2 Bit, 
which is equal to (2T) f 5) df, where A is the surface wave 
A 


fe) 
. ; : 2 
length. Integrating. Figure 3.10 yields © = .018. Using the 
We ie: _ 2 Bye ; 
approximation given in Beckmann, L” = 20° /L”, where o is obtained 


from integration and L is the correlation distance (3.2 cm for 


up-down wind and 10.2 cm cross-wind as shown later in this section), 


iz 


cross’ aad z* . Dd a 


2y were calcu- 
cross 


Z 
the values of % 5 ee 
ud 
lated and compared with the optically measured values of Ball and 
Carlson in Table 3.1. The relatively weaker cross-slopes in the 
present research is attributed to the completely covered water 
surface which maintained the wave build-up predominately in the 


windward direction. The reduced total slope is probably due to 


increased drag due to that greater covering. 


Table 3.1 


rms slopes of the surface (given by >= 20“ /L7) 


> 2 Z 2 
ud cross 
Scheible 
and .0096 .0008 .0104 
Fowler 
Ball and 
eageieos .0097 .0068 .01655 
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Analysis of the Energy Spectral Density 

It is believed that the surface height can be described as 
h(t) = E(t) Cos @(t). If this function could be determined, many 
powerful mathematical tools could be applied to thoroughly analyze 
the surface. That the surface can be described by h(t) is hinted 
by the fact that the energy in a given frequency band is time 
varying. Attempts were made to find E(t) and @(t) by studying the 
amplitude variation of the energy centered about a given frequency 
t and by studying the frequency modulation of the surface. 
Figure 3.11 shows the surface analog as measured at station 1 C. 
Et is obvious that there is some amplitude and phase/frequency 
variation. 

Amplitude Modulation, It is thought that the amplitude 
variations of a given frequency component of the surface (ft) 
may have a direct bearing on the scattered sound pressure. To 
find a description of these variations. Various experimental 
attacks were tried. First the wave anlyzer was set to 4 and 
the filtered output was analyzed. The three available filter 
settings (f. = 3, 10, 50 Hz) were used. Since the amplitude 
variations were at a very low frequency, the surface analog 
was replayed at 37.5 ips (tO x 1OU)> Figure 3,12 through 3.16 
show the time variations of given f,'s and Af 'S; in these 
figures all times and frequencies have been reduced to real values. 

The output from the wave analyzer was processed as shown in 
Figure 3.17. The output from the second tape recorder was cali- 


brated to the output of the first tape recorder (surface analog). 
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Figure 3.13 Wave Analyzer Output (f9=3.8H2,Afo=.lHz 
.001 sec/cm) 





Figure 3.14 Wave Analyzer Output 
(fo=3.8Hz,Afo=.5Hz, .001 sec/cm) 
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Figure 3.15 ave Anal . out 
(f.=h.8Hz, fy=.5 ,001 sec/em) 





Figure 3.16 Wave Analyzer Cutput 
(fo=2.8Hz, fp=.5Hz, ,OO01 sec/em) 
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Figure 3.17 Circuit to Measure the Energy Spectral Density 
of a Filtered Surface Component 






4f i= 3 Hz 


Figure 3.18 Circuit to Measure the Wave Analyzer Filter 
Response 
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The second wave analyzer was used with At, = ,.003 Hz in real fre- 


quency. The calibration was that the energy spectral density out 
of the second wave analyzer was equal to (mv) * x07 Xx 10° em? eee 
referenced to 1 Hz Bw. 

The theory behind this analysis is as follows. If a pure 
tone is amplitude modulated, the modulation spectrum will appear 
symetrically on either side of the tone in the frequency domain. 
Therefore, if one filters fo and measures the spectrum in the 
vicinity mn fo? the amplitude modulation will appear as side bands. 
This is exactly what was done; the results appear in Figures 3.19 
through 3.22. 

To verify that these results were not just functions of the 
filter in use, filter response was measured using the circuit in 
Figure 3.18. The SG-85 generates a carrier frequency t at 10 kHz, 
210 mv. This carrier is externally modulated by the HP 3300A. 
The output from the SG-85 is an amplitude modulated signal. The 
modulation frequency (f,) was varied in steps and the output from 
the SG-85 was fed to the first wave analyzer, set to 10 kHz, 
Af, = 50 Hz. The second wave analyzer was varied to measure the 
output of the first. By holding the amplitudes of f and f, con- 
stant, it was possible to determine the filter characteristics of 
the first analyzer. The result is shown in Figures 3.23 and super- 
imposed on the Figures 3.19 through 3.22. 

Aithough the theory is correct, it was not possible to find 


the spectrum of the on component, because f is not a pure tone 


but a component in a band of frequencies. It is not possible for 
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the filter to pass fo and its amplitude modulation alone; it must 
pass the components of the surface near fo: This is shown by the 
fact that the results in Figures 3.19 through 3.22 are nearly re- 
productions of the energy spectral density near fo: 

The frequencies near to that the filter passes tend to mask 
the amplitude modulation of f° As the bandwidth is increased, 
the amplitude modulation appears to increase (Figures 3.20 and 3.21), 
a result of the added frequencies being passed. To gain additional 
insight into the time-varying surface characteristics a digital 
analysis technique, with its inherent flexibility of filter speci- 
fications, was applied. 

To find the description of the amplitude variations of speci- 
fic frequency components of the surface, the surface was sampled 
2500 times at At = 10 ms intervals using the CDC 160 computer (2500 
samples is the approximate limit of the CDC 160 buffer). This block 
of 2500 was subdivided into blocks of m samples. Each block was 
Fourier analyzed using the Harm routine on the IBM 360 computer. 

The Harm routine determines the energy of specific components of 
the signal which are given by nAf, where n= 1, 2, 3, **** and 
Af = 1/mAt, where m must equal oe Tie highest frequency com- 


ponent that can be analyzed is given by BW = 1/2At. 


Table 3.2 


Studied Sample Blocks for Harm Routine 


At m Af BW 
Case l 10 ms 32 3.12 Hz 50 Hz 
Case 2 10 ms 64 1.56 Hz 50 Hz 
Case 3 10 ms 1024 .0975 Hz 50 Hz 
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The surface was examined for the three cases listed in Table 
3.1. The spectra obtained for the first two cases were not useful 
because they were not well enough defined over the region of in- 
terest. We are interested in the behavior of the surface between 
1 and 5 Hz. Case 1 gave one data point in this range, and Case 2 
gave 3. Case 3 gave many data points between 1 and 5 Hz. 

Because Case 3 required 1024 samples, only 2 separate spectra 
could be obtained from the 2500 samples. Comparison of these two 
spectra indicated that the energy was not constant in any given 
component. Whether this energy variation was a result of second 
order interactions between the components, a statistical variation, 
or a combination of both was not determined, 

It was hoped that one could select a frequency component and 
analyze its amplitude variations by taking its value in each of 
the consecutive subdivisions and performing:a second order Fourier 
analysis, where AT = mAt, BW = 1/2AT. This procedure was not con- 
sidered valid for cases 1 and 2 because of the lack of definition 
in the first order spectra. Case 3 was not valid because only 2 
points were available. 

It iS proposed that one could select At = 100 ms, Af = .1 Hz, 
BW = 12.8 Hz, and m = 256 and determine 10 consecutive spectra 
from 2560 samples. This may lead to greater insight into just how 
the energy varies in and between specific components. 

The variance of the wave analyzer output was determined for 
= 3.8 Hz and Af. = .03, «1, wo, > Hz using the circuit 
shown in Figure 3.24. The 5 Hz bandwidth was obtained by replaying 


the surface analog (tape recorded #1) at 3.75 ips (tf. x 10) and 


yl 


eetee f xl000 
T.Re-l W.A. T.R.-2 DETECTOR 


Figure 3.2 Circuit to Measure the Variance of the Output of the 
Wave Analyzer 
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setting Af, = 50 Hz. The analyzer output was recorded on the 
second tape recorder. The gains of this recorder and the wave 
analyzer were set to one. This set the output from TR-2 equal 

to the input from TR-1. TR-2 was replayed at increased speeds 
such that the output was fy x 1000. This output was fed to the 
Ball and Carlson sampler. The samples were recorded on the pulse 
height analyzer. These samples (Figures 3.26 through 3.28) were 


processed in the IBM 360 to obtain the variance, 


The results are plotted in Figure 3.25. 

Figure 3.25 shows the behavior of the filter. As the BW 
is decreased, the variance approaches zero; as the BW is increased 
the variance approaches the value one would have obtained if the 
wave analyzer had not been used and the surface had been processed 
directly. Notice, this would not be a measure of the rms height 
of the surface, but a measure of how the rms height changes. 

The PHA data, shown on Figures 3.26 through 3.28, indicate 
that the amplitude of the 3.8 Hz surface component has a stati- 


-0O3 Hz. 


iI 


stical distribution which is close to Rayleigh for Aft, 
As the bandwidth is increased the distribution changes from 
Rayleigh and becomes close to Gaussian for At = "soe 

Frequency Modulation. The surface analog oor x LOO) was fed 
into the HP 5210-A frequency meter, Figure 3.33. This meter in- 
dicates the average frequency of the signal. The needle on the 
meter varied from 3.2 to 4.8 Hz (real frequency) in an apparently 
random manner, indicating the surface is in fact frequency 


modulated, 
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Figure 3.27 PHA Data for Variance Measurements 
( f, = 3.8 Hz, Of, = ol Hz ) 
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Figure 3.28 PHA Data for Variance Measurements 
( fF, = 3-8 Hz, Af, = .5 Hz) 
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The frequency meter has an internal discriminator network 
which gives a DC output voltage proportional to the instantaneous 
frequency. This permits measurement of frequency modulation com- 
ponents up to 25% of the carrier frequency. For proper separation 
of the modulation from the carrier at the output, an internal low 
pass filter must be set less than, or equal to, 25% of the carrier. 
In this case a cutoff frequency of 1 Hz (real frequency) was used, 
.25 x 4 Hz. Figure 3.26 shows the input and output of the fre- 
quency meter. In this figure the lower trace is the instantaneous 
variation of the input frequency from the carrier. Figure 3.30 
shows a case where the input frequency does not vary in time. 
Figures 3.29 through 3.32 were obtained when the surface was 
filtered through the wave analyzer (f = 3.& Hz, Af, = -03, «1, 

.5 Hz) before going to the frequency meter. 

The discriminator output, representing the frequency modu- 
lation of the surface, was analyzed using the wave analyzer, 
Figure 3.33. The results are plotted in Figure 3.34. This is the 
spectrum of the signal which frequency modulates the surface. The 
exact mechanism which causes this modulation 1S unknown. The only 
thing that can be said is that it is relatively strong for modu- 
lation frequencies below .25 Hz; that is, the frequency modulation 
is very small compared to the predominate frequency in the surface. 

Conclusions. The wind-agitated water surface is a very com- 
plicated, non-linear, dispersive medium. Results obtained from 
our analysis of the surface show that the energy in the surface 


randomly drifts through the spectrum. At one instant it is 
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Figure 3.29 Discriminator Input (upper) and Output (lower) 
for the Unfiltered Surface Analor 
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Figure 3.30 Discriminator Input (upper) and Output (lower) 
for Filtered Surface Analog 
( fo = 3.8 Hz, Af, = .03H2) 
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Figure 3.31 Discriminator Input (upper) and Output (lower) 
for filtered Surface Analog 
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Figure 3.32 Discriminator Input (upper) and Output (lower) 
for Filtered Surface Analog 
( f, © 3.8 Hayat, = 5 Hz) 
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Figure 3.33 Circuit to Measure Frequency Modulation of 
the Surface 
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Figure 3.34 Surface Frequency Modulation 
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concentrated at 2.6 Hz) latermatmoeoeH: .fthen at 3.5 Hz, and so 
on. As this energy iS moving about, it is amplitude and phase 
modulating the surface components. 
Surface Wave Height and Distribution 
Electronics. Figure 3.35 shows the circuit used to measure 
the height distribution of the wave surface. All of the components 


of this circuit have been described earlier except the sampling 





circuit and the nuclear pulse height analyzer. Both of these were 
originally used by Ball and Carlson in their measurements of the 
surface distribution. 

Calibration. To calibrate the circuit so that the channel 
numbers on the pulse height analyzer could be directly related to 
wave height, the depth of the probe was varied with the water sur- 
face smooth. For each depth the output was sampled and read on 
the pulse height analyzer. In this way the entire circuit was 
calibrated, and it was not necessary to calibrate any individual 
piece of electronic equipment. A sample calibration curve is 
given in Figure 3.36. 

Results. The wave height distribution was determined for all 
nine stations shown in Figure 3.2. The distribution at station 
1 C is given in Figure 3.37. The curves at the other stations 
were similar. The solid curve is a Gaussian curve with zero mean 
and the same variance as our data. Table 3.2 gives the values of 
rms height measured at each of the nine stations. The mean value 
of rms height was found to be 0.21 cm with a standard deviation 
of .04 cm. Because of this large variation, and because the same 


20% variation was noticed when readings were taken at one point at 
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Figure 3.36 Calibration of Surface Probe to PHA 
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different times, it was considered necessary to record the surface 
distribution and compute the rms height each time a measurement of 
the scattered pressure distribution was taken. In this way it was 
possible to ensure that the theoretical value with which we com- 
pared our data was computed using nS Sim tensors rms height. In 
addition to determining the variance, we measured the skewness and 
the kurtosis or "peakedness"., Skewness is defined as one-half the 
third moment divided by the third power of the standard deviation. 
And kurtosis is defined as the fourth moment divided by the fourth 
power of the standard deviation minus three divided by two. The 
skewness is less. than 0.1 and the kurtosis less than 0.4 for aim 
surface height distributions. For a true Gaussian distribution 


these would be zero. 


Table 3.3 


rms wave heights 


Station rms wave height in cm 
le O.197 
1R 0.161 
LL 05200 
2C OgzZ25 
2R Olas. 
AG 0.263 
3C Oe 20 
3R 0.183 
3L 0.266 
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Surface Correlation Function 
The surface spatial correlation function is defined as 
1 T 
R(d) -iim = if h(t,0)h(t,d)dt 


Tro 


al 


where h(t,0) is the surface height at a reference point, and 
h(t,d) is the surface height at a second point d cm from the first. 
To find this function, two probes ere placed in the surface, 
Their separation was varied using the apparatus shown in Figure 3.38. 
For each separation, R(d) was computed using the circuit shown in 
Figure 3.39. Oscillator frequencies of 10 and 30 kHz were used 
to reduce cross-talk between the two probes. The surface analog 
from each probe was processed in the same manner as that used to 
measure the surface spectrum. However, the DC level had to be 
more completely removed by the differential amplifier. A .2 volt 
DC component in either signal would introduce a 4% error in the 
correlation function. Because only one differential amplifier was 
available with an internal low pass filter, the tape recorder was 
used to filter the 10 and 30 kHz ripples from the surface analogs. 
The tape recorder operating in the FM mode has an inherent high 
quality low pass filter which can be set to 100 Hz, 1 kHz, or 
10 kHz, depending upon the speed in use. The analog computer per- 


formed the following integrations: 


a 2 is 2 
Pao 0h(t,0)] dt Bei tad) ) dt 
it 
c =f n(t,0)h(t,d)at 
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where A = auto-correlation from probe l 


OS 
H 


auto-correlation from probe 2 

C = cross-correlation between the probes 
A, B, and C were voltages proportional to the integrations. The 
absolute calibration was not important because the cross corre- 
lation function was normalized by dividing by the geometric mean 
of the auto-correlations. Integration times of 30 seconds, one 
minute, three minutes, and five minutes were tried. The normalized 
correlations function was consistent for all times. However, the 
auto-correlations displayed slow temporal and spatial variations 


of as much as 20%. It is believed that these variations are due 





to second-order interactions between waves. 

Before data were recorded at each separation, the two auto- 
correlation were simultaneously obtained. These values should 
be the same. A substantial difference indicated either a DC 
component was present in one or both of the signals, or the gains 
were different between the two circuits. These two parameters 
were adjusted until the auto-correlations were within 5% of each 
other. An error of 5% in the auto-correlations would be reflected 
as an error of 5% in the cross-correlation for d less than 1 cm. 
As d is increased, C decreases because of the reduction of cross- 
correlation. However, A and B are not changed. Therefore, any 
error in A or B would cause an increased error in the normalized 
cross-correlation. This erroris variable, changing as the sepa- 
ration changes. Therefore, an error of + .1 in correlation was 


assumed for separations greater than 5 cm. 
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The correlation function was determined in the up/down-wind 
direction and the cross-wind direction, They are shown in Figures 
3.40 and 3.41. The correlation distances for these two directions 
are 3.2 cm and 10.2 cm, respectively. These results indicate that 
very little surface energy was propagated perpendicular to the 
wind direction, 1.e., ier K. 


y 


The theoretical developments consider the normalized surface 


‘ : ; -d ia Z 
correlation function p to go as either e /L or e (d/L) 


where L 
is the correlation distance. The up/down-wind function, at first 
glance, suggests a (Sin x)/x behavior. To show that it wasn't any 


ee 
"= a ee) , ana | (Sin 2e2ayar2|7 


of these simple relations, 
(2.2da/3.2) were plotted on Figure 3.37 for comparisons. 
The rms slope has been shown by Beckmann to be ) = \f20 /L. 


For the up/down-wind slope this would give an rms slope of ap- 


proximately ae 
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a The Specularly Scattered Sound 
General 

In order to measure the effects of the surface on sound wave 
propagation, sound was directed at the surface from one transducer 
and received by another. These transducers were mounted on lo- 
cally manufactured racks which have externally adjustable angle 
dials so that the angles of incidence and reflection can be ad- 
justed to an accuracy of better than ae 

Spectrum 

Underwater sound at 30 kHz was specularly scattered from the 
surface described in section 3. The angles of scatter were 9 = 21, 
35.3 and 51 degrees. The scattered pressure was observed to be 
amplitude modulated at a low frequency, on the order of 3 Hz. 
This signal could not be directly analyzed because the narrowest 
filter on the wave analyzer was 3 Hz. At 30 kHz center frequency, 
the spectrum appeared as a pure tone; that 1s, the modulation fre- 
quencies were masked by the carrier. If the frequencies could be 
multiplied before analysis, thus spreading the frequencies, the 
spectrum could be measured directly, as was done in the surface 
analysis. However, the wave analyzer had a maximum frequency 
response of 50 kHz. To circumvent this problem, the envelope of 
the scattered sound pressure was detected. This demodulated signal 
was analyzed using the same techniques used for the surface ana- 
lysis, Figure 4.1. The scattered sound pressure was also analyzed 


for frequency modulation. 
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Ampiitude Moduiation. The spectral density of the amplitude 
modulation was obtained for 9 = 2 ichae This spectrum is shown in 
Figure 4.3. In this case, the mean scattered pressure was 66 
microbars, the surface roughness was g = .16, and the energy 
density is expressed in Wea ae This absolute analysis 
was done for 9 = 2634 only. The spectra at 9 = 21° and 5 are 
in relative units; an absolute calibration was not performed. The 
spectrum at 35 degrees was thoroughly analyzed to find its re- 
lationship to the surface spectrum; the spectra at 21 and 51 degrees 
were obtained to confirm the findings. 

Frequency Modulation. The circuit shown in Figure 4.5 was 
used to measure the frequency modulation of the scattered sound 
pressure. To ensure accuracy of the experiment, a frequency syn- 
thesizer was used at a frequency of 30,000.0 + .1 Hz. To analyze 
the received signal, the HP 5210 frequency meter, which was des- 
cribed in section 3, was used. The cutoff frequency of the output 
filter (1 kHz in this case) must be chosen such that the carrier 
does not introduce ripple in the output. When a pure tone is 
put into the meter, the output is a DC voltage between O and 1 volt. 
The actual level depends on the frequency and meter range in use. 
When the carrier is frequency modulated, the output varies about 
this reference DC. The amount of variation is proportional to the 
frequency deviation, and the frequency of the variation is equal 
to the modulating frequency. 

In order to calibrate the system, sound pressure was specularly 
scattered from a smooth surface (9 = 35 0m The synthesizer fre- 


quency was changed + 10 Hz and the output ‘DC voltage was found to 
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deviate + 10 mv. This procedure was repeated for several frequency 
changes down to + .5 Hz. When the frequency change was + 1 Hz, 

a deviation of 1 mv was detected in the output. However, for 
frequency changes below 1 Hz no deviation was detectable above the 
noise. These deviations were in agreement with values calculated 
following the instructions in the HP 5210 manual. These results 
indicated the frequency deviation in the analyzed signal must be 
greater than 1 Hz to be detected. 

The surface was roughened and the sound scattered from it. 

It was impossible to detect any change in the output from the fre- 
quency meter. This result led to the following conclusions. 1) If 
there was a frequency shift in the carrier it was less than 1 Hz. 

2) If there was any frequency modulation the magnitude of the de- 

viations was less than 1 Hz. 

Analysis. Examination of Figure 4.3 leads to the following. 
At 9 = 35.3 degrees the spectrum of the amplitude modulation had 
a form which closely resembled the surface spectrum, except for 
the dips near the peak and the added peak at about 6 Hz. The 
problem was to find the reasons for this behavior. The first 
thought was doppler shifting of the pressure upon scattering, al- 
though Parkins and Wetzel had shown that doppler shifting should 
not occur in the specular direction. 

It was discovered that this spectrum was in fact a result of 
doppler shifting. Figure 4.6 depicts the scattering geometry in 
the vertical plane of incidence defined by the axes of the two trans- 
ducers and the perpendicular to the surface at the center of the 


scattering area. Consider an object as shown. When sound pressure 


on 


Arjouloey Butzeiyeosg «=g* ty amaty 
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is incident on this object, the majority of the energy will be 
scattered along path B. However, there will be diffraction side 
lobes to this scattered pressure. 

Consider the case in which the object moves with velocity V 
and, at a given instant is positioned such that the side lobe is 
coincident with the receiving transducer. The pressure in this 
lobe can be shown to be doppler shifted by the scattering object. 


From geometry the amount of doppler shift is given by 


f,=7 [Sin(9® - B)- Sin (9 +a@)] 


IK 


Considering @ - B and both small, this equation reduces to 


IR 


f + (B - @)Cos @ 


28 Cos 9 


rh 
\ 


mete, £Or specular scattering, —p = O and £. = 0; that is, this 


d 
agrees with Parkins and Wetzel. 
Wetzel and Parkins, considering a resonant interaction between 


the acoustic wave length A and the surface wave length A, show the 


diffraction grating equation is 


where p is the order of the lobe, 9, is the angle of incidence and 


Ji 
o,, is the angle of reflection. Letting 8, = 0 - 8, and 0, =O+daQ, 
and considering p = 1; that is the first side lobe, this equation 
reduces to ‘ 


~ 28 Cos 0 
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Ty a= fort is considered the velocity of the diffraction grating, 
that is, the velocity of the surface wave component having a wave 
length A, 

A 


— - ee ‘wena Cos e 


which becomes f, = f..., when A = /2B8 Cos 0 is inserted. 
This development shows that the doppler frequency shift is 
equal to the frequency of the surface component causing the shift. 
To have a shifted component present in the received sound certain 
conditions must be met. There must be a path BB! which satisfies 
the diffraction equation in such a way that the first side lobe 
falls on the receiver. This implies that for each doppler shifted 
component there is an appropriate corresponding angle £ in the plane 
of incidence (and similar angles outside this plane). No more than 
one shifted component can arrive along a given BB! path. To change 
fo: ie ¢ must change; then the relationship A = 4/28 Cos 9 no longer 
holds; and the side lobe moves off the receiver. 
When two tones arrive at the transducer slightly different in 
frequency, they combine to produce amplitude modulation at the 
beat frequency. The frequency of the combined signal is somewhere 
between the two tones, depending on their relative strength. The 
model being considered deals with the spectrum of the amplitude 
modulation, that is, the spectrum of the beat frequencies. A given 
beat frequency equals the amount of doppler shift fae 
As reported, it was not possible to determine directly the 
frequency spreading of the carrier which resulted from this inter- 


action; the measurement of the amplitude modulation spectrum was 


conducted in place of the direct measurement. 
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The dips near the peak of the spectrum in Figure 4.8 are 
difficult to explain. They seem to be a result of interference. 
For the geometry used in this experiment, these dips occur when 
there is a (n + 4)k, n= 1, 2, *°* path length difference between 
paths BB' and AA. For f, = 2.6 Hz this difference is 4A; for 
fi = 3.2 Hz the difference is 3/2\(A = 5 cm). This relationship 
would explain destructive interference if the signals arriving 
along the two paths were coherent. However, it has been shown 
that they are incoherent. There is a possibility that there are 
surface elements outside of the plane of incidence which satisfy 
both the angle relation (for the proper frequency) and the path 
difference requirement for interference. The geometrical relations 
necessary to test this speculation have not been developed. 
Another possible explanation for the peaks and troughs is that 
they are due to higher diffraction orders. 

The secondary peak that occurs at f, = 6°Hz, Fidume 438 is 
attributed to the effect of the transducer side lobe, This is 
discussed under results and conclusions. 

‘Study of 

Results and Conclusions. AThe surface energy spectral density 
E(f_); the scattered sound amplitude modulation energy spectral 
density C(f 4); the sound pressure level at the surface where the 
Path BBiescatters, SPL(B), shows that the shape of E(f,) 
resembles the shape of E(f\); but is displaced; the peak of €(f5) 
moves to the left (lower frequency) as 9 increases, as predicted 
by the doppler shift equation ti ~ 2B Cos 9; for a given fre- 
quency, the SPL(fB) becomes smaller as 9 increases, again as pre- 


dicted by the doppler shift equation; and the secondary peaks in 
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E(f,) at 6 Hz, 9 = 21 and 35.3 degrees, correspond to doppler 
shifted energy arriving at an angle which corresponds to the side 
lobe of the transducer radiation pattern. Not shown on these 
figures, but found from the sound pressure magnitude studies, is 
the fact that as 9 increases up to about 75 degrees, the variance 
of the amplitude modulation not only moves to lower frequencies, 
but also decreases in magnitude as QO increases, 

An explanation of how the doppler shifted energy is produced 
and arrives at the transducer was presented in the analysis section; 
however, no attempt was made to explain the variation of the mag- 
nitude of the energy. It is contended that the energy in a given 
doppler shifted component is a function of both the energy con- 
tained in the surface component that causes the shift and the 
sound pressure level incident on the surface at the point where 


the shift originates. 


€(f,) = F E(£.) 


where F = FL SPL(8),A] 
A is the surface area that satisfies the diffraction condition. 

It must be emphasized that the explanation of the amplitude 
modulation spectrum has been based on the simplifying assumptions 
that (a) scattering comes principally from the water surface close 
to the plane of incidence and (b) the diffraction grating equation 
can be applied to the very limited surface area where the rays are 


incident at angle B + A where A << 8 defines a bundle of rays of 


common angle of incidence. 
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After completion of this research, a recent paper by Marsh 
and Mellen (1966) was discovered in which the amplitude modulation 
in specularly scattered sound was ascribed (without proof or 
detailed data) to interaction between the unshifted and the doppler 
shifted energy. 

Distribution of Scattered Sound Pressure 

Underwater sound at frequencies between 20 kHz and 100 kHz 
was specularly scattered from the model sea surface described in 
section 3 at angles of incidence between 35 and 85 degrees, with 
particular emphasis on angles greater than 70 degrees. 

The circuit diagram for measuring the distribution of the 
scattered sound pressure is shown in Figure 4.10. This circuit 
yields a distribution from the pulse height analyzer. When the 
calibration is applied, the distribution of sound pressure is 
obtained. 

Calibration. The calibration is done in still water. [In 
order to obtain the proper calibration, the output voltage of the 
receiving transducer is amplified a fixed amount (usually X2 or 
X5) by an HP 467A amplifier. The output of this amplifier is 
then measured uSing a vacuum tube voltmeter. The signal is sampled 
and measured by the pulse height analyzer. This determines the 


mirror reflection value of mean square pressure. Then the selector 


switch on the HP 467A is changed to variable and the output changed, 


read on the voltmeter, and read on the PHA. This calibrates the 
PHA to the output of the amplifier. The amplifier is returned to 
its original fixed amplification setting, and the calibration is 


complete. The resulting calibration line is a plot of transducer 
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output times a multiplication factor (two or five) versus channel 
number on the PHA. A sample calibration curve is given in Figure 
4,11. To convert this calibration from output voltage to sound 
pressure, the response curve for the transducer was used. This 
curve gives the output in db re 1 volt/microbar as a function of 
frequency. 

Measurements. With the surface roughened as described in 
section 3, the received sound signal was sampled at a sampling 
rate of one thousand samples per second, and the samples were 
recorded on the PHA for seven minutes. [In all, seventy-four sets 
of data were recorded and are provided in Table 4.1. For g less 
than 0.3 the theoretical value of relative mean Square pressure 
is e *% and it was this theoretical value with which we compar ed 
our data. Except for four sets of data, all of our g's were less 
than 2OR2 3. 

As mentioned in section 3, during the same time that sound 
pressure data were taken, the wave height analog was recorded. 

This signal was processed later to determine the rms wave height 
which existed during the period of the sound scattering experiment. 
A program (Appendix A) was used in the IBM 360 computer to 

provide, given the PHA distributions for the wave surface analog 
and the sound pressure analog and the necessary calibration factors, 
the proper values of the rms wave height, the skewness and kurtosis 
of the wave height distribution, the mean scattered pressure, the 
variance, skewness, and kurtosis of the sound pressure distribution, 


the roughness factor g, the theoretical value of relative mean 
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PHA Channel Number 
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Ficure 4.ll Sample Calibration Curve 
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square pressure e ae and the ratio of measured relative mean 


* 


Square pressure to ee, Figure 4.12 shows an example of the 


distribution of the specularly scattered sound pressure. 

Results. Figure 4.13 is a graph showing the ratio of the 
ene ee relative mean square pressure to the theoretical value 
e 2 as a function of the angle of incidence for various frequencies 
for g< 0.2. The dashed line at 1 corresponds to the theoretical 
value e %. The solid line is the theoretical value including the 
shadowing function of Wagner. 

Figure 4.14 shows the variance of the sound pressure as a 
function of the cosine of the angle of incidence for various 
frequencies. 

Figure 4.16 shows the variance at a particular anole sasme 
function of frequency. 

Conclusions. Figure 4.13 clearly shows that for the largest 
angles of incidence (smallest grazing angles) the measured values 
of relative mean square pressure are lower than the value predicted 
by the theory of Beckmann and higher than predicted using the 
shadowing function of Wagner. It seems clear that shadowing does 
occur; however, the calculation of Wagner overcorrects for the 
shadowing by ignoring multiple scattering and the wave nature of 
sound. Wagner's theory uses ray geometry and assumes that any 
“particular point on the surface is either illuminated or not 
illuminated. Because the wavelength of the sound is the same 
order of magnitude as the rms height of the surface, it should 
be expected that bending in the wave front will occur and areas 


which Wagner predicts will not be illuminated will in fact be 


partially illuminated. 
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Figure 4.15 Variance vs. Angle of Incidence for Various 


Frequencies 
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Figure ).16 Variance vs. Frequency for 72° Angle of 
Incidence 
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The measured variances were all equal to or greater than the 
variances predicted by Beckmann. This 1s to be expeeted. if 
independent random variables are added, the variance of the sum 
is the sum of the variances. Therefore, it is believed that the 
variance predicted by Beckmann is the minimum variance to be 
expected with any additional variances due to the electronics, 
shadowing, and multiple scattering adding to this predicted value. 
The dependence of the variance on aes © is seen in Figure 4.14, 
for angles less than 75 degrees. This dependence is in agreement 
with Beckmann and indicates that variances introduced by the ex- 
periment are nearly constant. 

Also it can be seen in the same figure and also in Figure 4.15 
that the variance increases rapidly for angles greater than 75 
degrees. It is believed that an additional variance due to the 
shadowing causes this behavior. As shown in section 3, the 
Gaussian distributed rms slope is approximately 5 degrees. This 
implies that 99% of the slopes will be less than three times this 
value or 15°. Therefore, when the angle of incidence is approxi- 
mately Wee (grazing angle of 15°), some of the surface should begin 
to be shadowed. The variance in the shadowing effect is then added 
to the other variances. 

Figure 4.16 shows the variation as a function of frequency 
for a particular angle of incidence aie It can be seen from 
this graph that the variance 1s approximately proportional to 
the second power of frequency. This is not in agreement with the 
prediction of Beckmann that the proportionality should be with the 
sixth power of frequency. No explanation for this discrepancy has 


been found. 
108 


It is interesting to notice that the various theories of 
sound scattering always asSume a Gaussian distribution of surface 
heights and, for g << l, predict a Gaussian distribution for the 
scattered sound. In fact, neither the real sea nor our model 
shows an ideal Gaussian distribution; and, of course, the scat- 
tered pressure is not an ideal Gaussian either. A virtue of 
the model experiment is that we can measure the resulting sound 
scatter distribution for a height distribution that shows re- 


alistic non-zero skewness and peakedness. 
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> Suite 

Experiments were conducted to measure the aneaeteen and 
spectrum of sound which was specularly scattered from a wind gen- 
erated model sea surface. The energy spectral density of the 


-9.3 


sur faces peaksmatue3.S Hegudnoppanametiects £ with increased 


-5.) 


frequency near the peak, becoming f for frequencies greater 


than 5 Hz. This frequency behavior is in agreement with the 
Neumann-Pierson spectrum and Phillips't theory which predicts fame 
dependence. The rms height of the surface was about 0.21 cm. 
The wave height distribution was nearly Gaussian. The spatial 
correlation function of the surface was determined for both up/ 
down-wind and cross-wind directions. The correlation distances 
were 3.2 and 10.2 cm, respectively. From the correlation dis- 
tance and the rms height, the rms slopes were calculated for up/ 
down and cross-wind directions. These values were .0096 and 
.0008, respectively. These results confirmed the assumption that 
the surface propagation vectors were nearly parallel down-wind 
with very little divergence. 

The distribution of the specularly scattered sound pressure 
was measured for various angles of incidence from 35° to grazing 
and various frequencies between 20 and 100 kHz, with particular 
emphasis on angles near grazing. Comparison was made with theories 
that predict the values of the mean and variance of the specularly 
scattered sound distribution. The mean values are in agreement 
with theory except at very near grazing indidence where shadowing 
causes less specularly scattered energy than predicted for no 


shadowing and more than predicted for geometrical shadowing. The 


Lie 


measured variances were generally greater than predicted due to 
the additative experimental variances. The observed dependence 
of variance on the fourth power of the cosine of the angle of in- 
cidence was in good agreement with theory. However, the variance 
was proportional to the second power of frequency rather than the 
sixth power as predicted. 

The spectrum of the specular scattered sound contains a 
predominant component at the frequency of the radiated sound. 
In addition it contains doppler shifted components, although 
theory predicts that the spectrum of purely specularly scattered 
sound should not contain doppler shifted components. 

The doppler shifted energy arrives at the wceiver from 
directions other than the purely specular. This is a result 
of finite divergence of the incident sound. It was found that 
the doppler shifted spectrum is directly related to the frequency 
spectrum of the surface. The energy level of the doppler shifted 
components is a function of the surface energy spectral density 


and the incident sound pressure at the off-axis points. 
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